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ABSTRACT 


An investigation into the limitations on the enhanced field-emitted current density 
in a fast-pulsed (rise-time = ns), high voltage (> 10°V), 1-inch vacuum diode was 


conducted using a computer simulation based on the Fowler-Nordheim equation. 


Oscillations in the emitted current density (due to the change in the amount of 
Space charge within the gap) were found to quickly decay into a final steady-state for 
the voltages applied. Steady-state values for a wide variety of work functions, electric 
field enhancement factors (based on the theory that "whiskers" on the cathode surface 
experience varying degrees of enhancement), and applied potentials were compared 
to two benchmarks: the amount of current density required to explode a whisker in 
< 10 ns by joule heating (Je = 10° A/cm*); and the Child-Langmuir (C-L) space- 


charge-limited current density. 


Steady-state values were found to be less than Je. One model of the formation 
process of a plasma at the cathode surface requires that Je be met or exceeded by 
the steady-state value. Thus, such a model is not supported by this project's findings. 

The C-L limit is based on a thermionic-type emission process. As only pure field 
emission (i.e., no thermionic emission included) was considered, the steady-state 


values were, in all conclusive cases, less than the corresponding C-L limited values. 
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l. INTRODUCTION 


The breakdown of the surface of the cathode in a pulsed vacuum diode has been 
the subject of much investigation in recent years. The existence of crater-like cathode 
spots after breakdown can be explained by unipolar arcing, but the phenomena 


occurring prior to the arcing are still debated. 


One model contends that a metallic "whisker" on the surface of the cathode 
provides a source for the enhanced emission of electrons. After a period of time (less 
than 10 ns for a fast-pulsed diode), the whisker explodes (due to the rapid build up of 
energy produced by the emission process) into a dense plasma which in turn drives 


the unipolar arc. 


This project will attempt to show that field emission alone (i.e., not including 
thermionic emission) from a whisker cannot reach high enough values to trigger the 
explosion mechanism. A computer simulation utilizing the Fowler-Nordheim equation 
for field emission will be used. The results will also be compared to the well-known 


Child-Langmuir space-charge-limited current density for a vacuum diode. 


A look at the dimensionless variables used in the program and the program itself 


are included in the appendices. 


ll. BREAKDOWN AND EMISSION THEORY 


A. BACKGROUND 


For over six decades scientists have been investigating pulsed-voltage vacuum 
diode breakdown. With the advent of newer technologies, shorter and shorter 
voltage-pulse rise times have become possible, with values ranging from 10° seconds 
in 1928 to 10°° seconds today [Ref. 1]. Corresponding to today’s very short rise times 
has been the discovery that breakdown even occurs on the nanosecond time scale. 
The effects of the breakdown process on the surface of the cathode have been well 
documented [Refs. 1-3], and a theory of the formation of cathode spots (Fig. 1) via 
unipolar arcing’ has recently been proposed [Refs. 2,3]. In dispute, however, is the 


process by which the plasma driving the unipolar arc is formed. 


Cathode spots are characterized by the formation of a dense quasi-neutral 
plasma near the cathode surface. There are two basic schools of thought concerning 
the primary mechanism for breakdown and the source of the plasma, only one of 
which will be attended to here. Anode-initiated breakdown would occur when ions 
released from the anode (due to bombardment of the anode by electrons emitted from 
the cathode) interact with the cathode to form the plasma. Here _ fast-pulse 
( = nanoseconds) breakdown shall be investigated, in which the anode-emitted ions 
would not have time to reach the cathode (for the 1 in = 2.54 cm diode gap of interest) 
before breakdown occurs; anode-initiated breakdown shall be discussed no further.2 


' The concept of an arc which originates and terminates on the same surface (i.e., “unipolar") 
was first expounded by Robson and Thonemann in 1958, and was later expanded by Schwirzke 
in the 1980’s [Refs. 2,4]. 


2 For a more thorough treatment of anode-initiated breakdown, see [Ref. 1]. 


Figure. 1. 


I 
® 100 ym 


Crater-like cathode spots formed by unipolar arcing. 


In the cathode-initiated breakdown school, it is proposed that all of the ions and 
electrons which make up the plasma are produced by processes which are initiated by 
events occurring only at the cathode [Ref. 1]. Within this school, however, there are 
many camps, each with its own views on what processes do indeed occur, and which 
events really take place to initiate them. One of the more reknown camps is led by 
Schwirzke, who recently proposed that neutrals released from the surface of the 
cathode are ionized by the emitted electrons, thus forming a plasma near the cathode 
surface (Refs. 2,3]. Another is led by Dyke (et al)? who proposed that joule heating of 
an electron emission site causes the site to vaporize, and in the process ions are 
created which allow the plasma to form [Ref. 1]. Here shall be examined a third camp, 
championed by Mesyats [Ref. 1], in which Dyke’s theory is taken to its most extreme 
limit, that of an explosion of the emission site into a plasma directly. The goal is to 
show that an elementary version of a Mesyats-type "explosive" model cannot support 


breakdown on a time-frame indicative of fast-pulsed voltage vacuum diode breakdown. 


B. THE MODEL DIODE 


The model's anode and cathode will consist of two polished metal plates 
separated by a vacuum. Despite the polishing and the degassing effect of the 
vacuum, however, some small-scale "irregularities" will nevertheless exist on the 
surfaces of the plates. Weakly bound bits of dust, oxides, and non-metallic inclusions 
litter the cathode surface, providing excellent sites for the enhanced field emission of 
electrons [Ref. 2]. These so-called "whiskers" will vary in size and shape, with the 


more needle-like (height-to-base ratio >> 1) metallic protrusions providing the most 


3 See [Ref. 1] for a list of references on Dyke's experiments. 


enhancement of the field emission. Flatter, non-metallic (e.g., dielectrics, semi- 


conductors, etc.) sites can also provide significant enhancement [Ref. 5]. 


C. ELECTRON EMISSION FROM THE CATHODE 


There are two mechanisms for the electric-field-enhanced emission of electrons 
from a metallic surface: thermionic (high temperature, low applied field); and field (low 
temperature, high applied field) [Ref. 4]. The presence of the electric field at the 
surface lowers the potential barrier at the plate-vacuum interface, thus allowing 
electrons in the conduction band an opportunity to become free from the metal (Fig. 2). 
In thermionic emission, electrons at high temperatures gain enough kinetic energy to 
pass over the barrier and into the gap. In field emission, the electrons tunnel through 
the barrier instead. In general, the electric field at the emitter surface must be on the 
order of 10’ V/cm for field emission to take place [Ref. 6]. For a diode gap of 2.54 cm 
(= 1 in), an applied potential of 2.54x10’ V would be necessary to achieve field 
emission, but as noted previously, the whisker emission sites experience an enhanced 
applied electric field (Fig. 3). The enhancement factor, m, varies from 1 to several 
hundred, with needle-like whiskers having the largest values [Ref. 7]. For electric 
fields less than 10’ V/cm, thermionic emission will dominate for temperatures above 
approximately 1500 K, with little substantial emission at lower temperatures [Ref. 6]. 

For the case of the non-metallic inclusions, a combination of field and thermionic 
emission occurs. Electrons tunnel from the metallic cathode material into the 


dielectric, then are thermionically emitted into the gap (see Fig. 4) [Ref. 5]. 
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Figure 2. Lowering of potential barrier at metal-vacuum interface (Ref. 7]. 
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Figure 3. Enhancement of electric field at whisker tip [Ref. 7]. 
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Figure 4. Combination of field and thermionic emission for a dielectric [Ref. 5]. 


For this elementary model, only field emission from metallic whiskers shall be 
considered, with the following assumptions: 
1. The time-scale for the explosion of a whisker is only a_ few 
(< 10) nanoseconds; 
2. The applied voltage rises instantaneously to a constant value at t = 0; 


3. The gap is perfectly evacuated; no positive ions are present, so no heating of 
the site by ion bombardment is possible; 


4. The initial temperature of the cathode (=300K) is too low for any 
appreciable amount of thermionic emission to occur at first; and 


5. Any increase in temperature during the short time-frame prior to explosion 
shall be neglected (and thus so shall thermionic emission) in this elementary 
model. 


Pure field emission is governed by the Fowler-Nordheim equation [Ref. 6] 


F2 ( -6.83x107 ) f(y) (oe 
J(F) = ( 1.54x10™ ) — e i. (1) 
Since f(y) (the Nordheim elliptic function) is slowly varying, it can be approximated as 


a constant* 
fy) = 0.94 , 
which leads to (see Figs. 5a,b) 


3/2 
2  (-6.42x107 ) 2 
J(F) = (1.54x10° ) 7 e En (2) 


where 


J = emitted current density (A/cm? ); 
= -mE = enhanced electric field at the whisker tip (V/cm ); 
m = electric field enhancement factor; and 


o = work function of the cathode metal (eV). 


J is a function only of F for any given whisker; any negative (i.e., electron) space 
charge which exists in the gap will reduce the emitted current density. If the applied 
potential remains constant, it can be predicted that J will start at some initial value, Jo , 
and decrease as more and more space charge enters the gap. The electrons, once in 
the gap, will be accelerated across to the anode, eventually disappearing there. When 
the first emitted electrons reach the anode, the electric field at the cathode will 
increase due to the loss of space charge. The increased field will cause an increase 


4 For a more thorough treatment of f(y), see [Refs. 4,6]. 
5 Eq. (2) requires that F be defined as a positive quantity, hence the negative sign. 
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Figure 5a. Log plot of J vs. F (Eq. 2). 
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Figure 5b. Expanded log plot of J vs. F for various work functions (Eq. 2). 
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in J, but soon the increased addition of space charge will exceed the loss at the 
anode; the amount of space charge in the gap will then rise, the electric field at the 
cathode will decrease, and J will fall. This oscillatory action should diminish as the 
emitted current density at the cathode comes to equal that going into the anode. In 
other words, the current density becomes constant at some final space-charge- 


determined value, Jr. 


It should be noted that in reality, the applied potential is not constant from t = 0. 
Rather, it ramps up from zero to some final value. The ramping has the effect of 
causing the current density to rise in a similar fashion; oscillations such as those 


described above are not observed. 


D. NON-APPLICABILITY OF CHILD-LANGMUIR LAW 


It is somewhat important to note, at this point, that Je is not the well-known 


Child-Langmuir space-charge-limited current density (A/cm?) 
V2 a 


de 


2e 


Mi ; (3) 








_ 4 
Jo. = 9 EQ 
with 


Eo = permitivity of free space = 8.85x107'* C*/N-m? ; 
e = lelectron charge | = 1.60x107'9 C ; 

M = electron mass = 9.11x10~°! kg ; 

V = applied voltage (V) ; and 

d = diode gap length (cm) . 


Je should be, in fact, somewhat less than Jo, . This difference arises due to the fact 
that Jc, is not derived for the case of a pure field emission mechanism. Rather, there 


must be an ample supply of electrons (eg., in the case of electrons emitted from a 
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heated filament) always available for emission. In the derivation of Jc,, one of the 
boundary conditions is that the electric field at the cathode go to zero due to the 
presence of the space charge (see [Fefs. 4,6,7]). If this were to happen to a pure field 
emitter, then J would go to zero, and a steady state would not develop. Since the 
field emitter cannot emit at Jc, for any appreciable length of time without shutting itself 


off, its steady state emission value must be something less: 


Je < Jer ( pure field emission ) . (4) 


E. EXPLOSIVE MODEL CURRENT DENSITY REQUIREMENT 


It has been shown by Schwirzke [Ref. 2] that, for stainless steel, a current density 
of Je = 109 A/cm? would be required to explode a whisker on the cathode surface in a 
few (<10) ns by joule heating. Thus, a whisker would need to emit a Jay = Je_ for 
approximately 10 ns before exploding. If the transient oscillation period is small 


compared to 10 ns, then Jayg = Jf, and 
Je 2 Je (explosive emission) (5) 


is required for the explosive emission model. This project will attempt to show that, for 
a wide range of applied potentials, work functions, and enhancement factors, Eq. (5) 
cannot be satisfied, thus proving that field emission alone is not capable of providing 


the current density required by the explosive model. 
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Hl. SIMULATION 


A. BACKGROUND 


The type of simulation model used here is, as [Ref. 8] so aptly puts it, one of the 
"simplest" possible, viable since the first digital computers came into use in the 1950's. 
While almost primitive in comparison to the latest innovations in so-called "particle-in- 
cell" simulations [Ref. 8], the model used here nevertheless is time-tested and has 
worked quite well for the rather simple events under investigation. Although the basic 
techniques are over three and a half decades old, the actual program is original and is 
included (in one form) in Appendix B for scrutiny. The program is written in the "C” 
language and was run on a Sun SPARCstation 2 incorporated into a system running 


SunOS™ 4.1.1.. 


B. SIMULATION THEORY 


In general terms, the simulation works as follows. Space charge is simulated not 
as a continuous beam of electrons, but rather as a series of infinitely thin, numerous 
charged disks of varying charge density o [C/m?]. In each time interval, the forces on 
the disks within the gap are computed, the disks are moved, and a new disk is created 
at the cathode. When a disk reaches the anode, it is no longer used in any further 
calculations. In the final steady state, all of the disks will have the same charge 
density; for every time step just as much charge is lost at the anode as is created at 


the cathode. 
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1. Electric Field 


The charge density o is calculated by multiplying the time step At by the 


current density J (Eq. 2)® 


Co = =—Ja4l. 


(6) 


The field-emitted current density at the cathode surface depends on €,,(0), 


the value of the electric field at the cathode due to all of the space charge disks, and 


Esp» the component of the electric field at the cathode surface due to the applied 


potential across the diode. In order to find E,.(0), one starts with the one-dimensional 


Poisson's Equation [Ref. 6] 


dic V(x) = q(nj — Ne) 


dx? Ep 


which, for the ion density n; = 0 and the electron charge q = -e becomes 


d * V(x) _ @ n(x) | 
dx? Eo 


where 


V = potential (V) ; 

x = position (= 0 at cathode, = d at anode ) (m) ; 
e = lelectron charge! = 1.60x107'9C : 

n = electron number density (m~*) ; and 


Eo = permitivity of free space = 8.85x107'* C2/N-m? . 


5 The current density J will always be defined as positive, hence a negative sign has been ad- 
ded on the right side of Eq. (6). 
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(7) 


(8) 


Using dV(x)/dx = —E(x), one gets 


GE(x) _ -€ n(x) (9) 
dx Eo 


Integrating both sides over the gap, using the relation E,.(d)=—E,,(0), and 


rearranging yields 


{ d 
E..(0) = Beo Jo e n(x) dx . (10) 


Now, using the charged disk concept, the integral becomes a summation 


1 
E = — Ax. 11 
sc (0) De zen x (11) 
where __n, = electron number density [m~°*] in disk number k. Realizing __ that 
enAx=-o yields 
1 
E..(0) = — ; 12 
sc ( ) 2€5 Pp —O, ( ) 
Finally, using Eq. (6), 
1 
E = — At . 1 
sc (0) 2€, p Jk t ( 3) 


For a disk located at a point x’ in the gap, E,,(x’) can be determined from 


Eqs. (9) and (10) by breaking the integral of the left-hand-side into two parts and 


subtracting: 
Esc (X) ae 
J joe eo [, -@ nx) dx (14a) 
becomes 
Esc(X’) - Ese(0) = -> F dat, (14b) 
£0 0x’ 
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and 


Es.(d) , 3 
dE,.(x) = — |, -e n(x) dx (15a) 
ald ies Eo Je 
becomes 
Pad) s(x) Ate (15b) 
£0 xd 


Subtracting Eq. (15b) from Eq. (14b) and rearranging yields 


, 1 
E.. (x’) = De XY At = a > J, At ° (16) 
x’ 


aR 
Eq. (16) is equivalent to saying that the space charge "in front" (x’-d) will contribute 
positively to the electric field at x’, whereas the space charge "behind" (0-—x’) will 
contribute negatively. 

If the applied voltage remains constant, so will the applied electric field E.p5. 


The total electric field at x’ is therefore’ 


, 1 1 
ET(x) = —Eapp 3 2€ bk At om Deo », J, At. (17) 
x’ 


0x’ 


At the cathode, the enhancement factor must be included: 


pilliae 


F = —m E+(0) = mE = 2€, 


y at} . (18) 


gap 


A schematic diagram of the electric field components in the diode is shown in Fig. 6. 


7 Like F, Ego will be defined as a positive quantity, so the appropriate sign change must be 
made. 
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Figure 6. Schematic diagram of electric field components in diode. 


ty, 


2. Equation Of Motion 


The nature of the large (210° V) applied potentials used in the simulation 
requires that relativistic equations of motion be used. For any electron in a disk ata 


position x, the relation is (based on mass x acceleration = charge x E-field) 


ayB) = 52 Eyx) = <2 Entx) , (19) 
where 
v= 1/ a — Vv" / Co: 
B=Vi/c; 
t = time (s) ; 


e = lelectron charge | = 1.6x10°'9 C ; 

M = electron mass = 9.11x107! kg ; 

c = speed of light in vacuum = 3.0x10° m/s ; 

E;(x) = total electric field at position x (V/m) ; 

xX = position ( = 0 at cathode, = d at anode ) (m) ; and 


V = velocity (m/s) . 


Substituting into Eq. (19) and carrying out the derivative gives 


vot + = <2 Exx) . (20) 


Taking the derivative of y and substituting in yields 
2 dV —e 

— = — : 21 

(Bey) ee (21) 


With a little more substituting, Eq. (21) can be solved for dV / dt: 


aq -e 
at Me E+(x) . (22) 





Eq. (22) is the equation of motion used in the simulation. 
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IV. RESULTS 


A. OBJECTIVE 


Again, the goal of the simulations is to show that the value of Jr for a given set 
of conditions will not be greater than the 10° A/cm® required by the elementary 


explosive emission model (Eq. (5)). 


B. METHODS OF OBTAINING RESULTS 


For lesser values (< 10° A/cm?) of Jo (the value of the current density at at the 
cathode at t = 0; j.e., with the potential applied but no space charge in the gap yet) 
the program is easily capable of running through the transient period and into the final 
state (as in Fig. 7a) in a short (= computer hours ) period of time. But for larger 
values of Jo, the run time goes up considerably (= computer days ). The objective 
requires only that Je <Je is proven; it is possible to do this without running the 
program through to the final state by using three simple "tricks": 

1. Attempt to simulate the steady-state electric field at the cathode by artificially 
adding in, at t = 0, an estimate of the steady-state space charge. The 
artificial charge will be an immovable source located at the anode so that it 


will always be “in front" of the emitted charges; 


2. Remove equal parts of the artificial space charge during each time step of the 
simulation, until it is all removed; and 


3. Set the amount of artificial space charge removed each time step equal to the 
amount created at the cathode in the first time step. 
Although the resulting plots will certainly be non-physical (no such artificial space 
charge exists in reality), the graphical output can nonetheless be _ definitively 


interpreted. 
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During the early phase of the removal period, the graph of J vs. t will be a 
horizontal straight line, as in the final steady state (just as much space charge is 
removed each period as is put in, so J will not change) (see Figs. 7b,c,d). If the 
estimate of the steady-state space charge is exactly correct, the graph of J vs. t will 
continue to be straight, even after all of the artificial space charge is removed, 
because in this case the simulation will have started equivalent to the steady state. 
Thus, if the graph remains flat after all the artificial space charge is gone, the value of 
Je will be known exactly (see Fig. 7b). Guessing the exact value of the steady state 


space charge is, however, highly unlikely. 


If the estimate of the steady-state space charge is too high, the value of F will be 
smaller than in the actual steady-state case, and the value of J during the early phase 
of the removal period will be lower than the actual Je. In this case some of the 
artificial charge will still exist when the first emitted electrons reach the anode. At this 
point, the amount of charge lost each time period will increase, thus increasing F and 
J (see Fig. 7c). 

More useful, however, is the case where the estimate of the steady-state space 
charge is too low. The value of F will be larger than in the actual steady-state case, 
and J will be larger than the actual Je. In this case, the artificial charge will be used up 
before the first emitted electrons reach the anode. Once all of the artificial charge is 
gone, the value of F will drop as more charge enters the gap while none is removed 
(see Fig. 7d). It can be concluded that the early removal phase value of J is greater 
than the actual Je. An upper bound on Jr can thus be found by this "maximum value 
method", and if the calculated upper bound is less than Je, the objective will have 


been met. 
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Figure 7a. Log of the current density ratio vs. time. No artificial space charge is 
added, and the simulation is run until the final steady state is reached 


(2 ns). The steady-state ratio is J/Jy = 1.37x107. 
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gap = 2.54 cm 

m= 10 

E.pp = 3.94e+06 V/cm 
= 1.68e+05 A/cm 


Total Time Steps = 1000 
Artificial Space Charge = 12.3752 


Log J/Jp 
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Figure 7b. Log of the current density ratio vs. time. Enough artificial space charge 
was added to simulate the steady state of Fig. 7a exactly (12.3752 
dimensionless units). The value of f is the number of time steps required 
to use up the artificial charge. The value of old indicates the time step of 
creation for the charge disk closest to the anode. Iff>0, old > 1, and the 
plot is straight, then the line is at the equilibrium value. Note the run-time 


is 0.1 ns (1/20 of the full-run time-steps in Fig. 7a). 


22 


gap = 2.54 cm 
10 
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Figure 7c. Log of current density ratio vs. time. More than the equilibrium value of 
the space charge has been artificially inserted (14 dimensionless units). 
In this case, f = 0, old > 1 indicates that the artificial space charge has not 
been used up before the first disk reached the anode; the plot rises at the 
end as expected. Again, only 1/20 of the full-run time-steps is required 


to obtain a result (in this case, a lower bound on Je). 
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Log J/J 


-3.00 


Figure 7d. Log of current density ratio vs. time. Less than the equilibrium value of 
the space charge has been artificially inserted (11 dimensionless units). 
In this case, f > 0 and old > 1, but the plot drops from the early-phase flat 
portion. 


steady state value (-1.863) of Fig. 7a. An upper bound on the current 
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density has been determined in only 1/20 of the full-run time-steps. 
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Note that the early phase value (-1.588) is greater than the 

















C. SIMULATION RESULTS 


As can be seen in the last column of Table 1, under nearly all tested conditions, 
Je<Je= 10° A/cm?. In fact, in most cases, Jo was less than Je, indicating that only 
for high applied electric fields (> 3x10° V/cm) and large enhancement factors (>100) 
would there even be a possibility for the final value of the current density to exceed 
the explosive limit. 

For those cases where it was not shown that Je <Je (denoted by ¢ ), the 
excessive amount of run time (= computer days) to lower the maximum level was just 
too prohibitive. 

It was predicted earlier that Je would be less than the Child-Langmuir space 
charge limited current density Jc, (Eq. (4)). The program was run long enough (when 
possible) for the @ = 3.0 eV case to determine if the maximum of Je is indeed less than 
Jc. The results are summarized in Table Il. Similar results are expected for 
d = 4.0 and 5.0 eV as well, since the values of Jp (and logically the maximum Je also) 


are lower for these two cases than for the @ = 3.0 eV case. 
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TABLE |. SIMULATION RESULTSt 


E.., (Viem) | (eV m_ | Jo (A‘cem? Je (A/cm? 
3.94x 10° 5.0 1 | =0 =() 

3.94 10° 5.0 10° | =0 =a 9) 
3.94109 5.0 100 | 5.86 4.93 

3.94x 10° 5.0 300 | 9.92x10® < 7.84x10& * 
3.94 10° 4.0 1 | =0 =0 

3.94 x 10° 4.0 Vor cate = 

3.94 10° 4.0 100 | 1.30x10° 8.69x10! 
3.94 10° 4.0 300 | 6.98x10’ < 6.97107 * 
3.94x 10° 3.0 1 | =0 =) 
3.94x10° 3.0 10 | =0 = 

3.94 10° 3.0 100 | 1.68x10° 2.17x10? 
3.94x 10° 3.0 300 | 4.27x10° < 5.00x10° * 
3.94 10° 5.0 1 | =0 = 

3.94 x 10® 5.0 10 | 5.86 5.77 

3.94 x 10° 5.0 100 | 7.73x109 < 6.81107 * 
3.94 x 10® 5.0 300 | 2.34x10'' | < 4.24x107 * 
3.94 x 10° 4.0 Je) ee =i) 

3.94 10° 4.0 10 | 1.30x10° 4.84x102 
3.94 10° 4.0 100 | 1.62x10'° | <7.24x107 * 
3.94 10° 4.0 300 | 3.48x10'' | <7.39x10’ * 
3.94 x 10® 3.0 1 | =0 =) 

3.94 x 10° 3.0 10 | 1.68x10° 2.30x10° 
3.94 x 108 3.0 100 | 3.42x10'® | <8.00x10® * 
3.9410 3.0 300 | 5.41x10'' | <9.00x10° * 
3.94x 107 5.0 1 | 5.86 5.85 
3.9410’ 5.0 10 | 7.73x102 < 6.81x10’ * 
3.94x 107 5.0 100 | 3.98x10'* | <8.00x108 * 
3.94107 5.0 300 | 4.05x10' | <7.00x10°9# * 
3.94 x 107 4.0 1 | 1.30x10 1.07x10° 
3.94107 4.0 10 | 1.62x10'° | <7.24x107 * 
3.9410’ 4.0 100 | 5.25x10'* | <8.00x108 * 
3.9410’ 4.0 300 | 5.15x10'? | <7.00x109# * 
3.94107 3.0 1 | 1.68x10° 1.55x104 
3.9410’ 3.0 10 | 3.42x10'° | < 8.00x10° * 
3.9410’ 3.0 100 | 7.32x10'* | <8.01x10° * 
3.94x 107 3.0 300 | 6.97x10'% | <7.00x109 * 


t For a diode with plate spacing of 1 in = 2.54 cm. 
* Indicates maximum value method result. 
+ Indicates run-time-constrained result. 
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TABLE Il. 


SIMULATION RESULTS VS. CHILD-LANGMUIR LAWt 


» (eV E... (V/cm) | Jc, (A/cm?) | Jo (A/cm?) | maximum Je (A/em? 

3.0 1 3.94107 3.62x10° 1.68x10° 1.55x104 

3.0 10 | 3.94x10® 1.14x104 1.68x10° 2.30x10° 

3.0 10 | 3.94107 3.62x10° S42 <4101° 8.00x108 # 
3.0 100 | 3.94x10° 3.6210? 1.68x10° 2.17x10 

3.0 100 | 3.94x10® 1.14x104 3.42x10'° 8.00x106 # 
3.0 100 | 3.94x107 3.62x10° 7.32x10'? 8.01x108# 
3.0 300 | 3.94x10° 3.62x10? 4.27x108 5.00x10° * 
3.0 300 | 3.94x10® 1.14x104 5.41x10"" 9.00x10& # 
3.0 300 | 3.94x107 3.62x10° 6.97x10'9 7.00x10° # 


D. OTHER GRAPHICAL OUTPUT 


The plot of Fig. 7a appears to agree with the basic theory that the value of the 
emitted current density will oscillate at first, then settle down until a final steady state is 
reached. But earlier versions of the program encountered problems. In an effort to 
determine the causes of the difficulties, the program was altered to produce plots of 
other variables. Figs. 8-11 are some examples run for the same conditions. Although 
the applied electric field is rather high in these plots, the results are nonetheless 
representative of the output for the applicable range of electric fields. Short computer 


running time was the primary factor in the choice of parameters for these examples. 


ft For a diode with plate spacing of 1 in = 2.54 cm. 
+ Indicates run-time-constrained resuk. 
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Fig. 8a shows how the velocity of the disks varies with distance. This plot shows 
that the need for relativistic equations of motion is quite real. Fig. 8b is a blow-up of 
Fig. 8a which shows that the velocity reaches 2.98x10'° cm/s by the time a disk gets 


only 5 % of the way across the gap. 


Fig. 9a shows how the density, similar to the velocity, is nearly constant across 
the gap (at equilibrium). This is to be expected at equilibrium, as all of the disks have 
the same charge density and most are moving at essentially the same speed. Fig. 9b 
is a blow-up of Fig. 9a, showing that the value of the density at 5 % of the way across 
the gap is within approximately 0.5 % of the value at the anode. The density is 
greater near the cathode because the slower moving electrons there are still 
"punched", whereas those farther away are fully separated by the acceleration 


process. 


If the density is essentially constant, the electric field should appear as a sloping 
Straight line centered about E,,,. Fig. 10a shows the magnitude of the electric field 
vs. distance. As predicted, the plot is straight and centered about E,,,. In order to 
see the effect of the higher density near the cathode, the plot was blown-up into 


Fig. 10b. Indeed, a small deviation can be seen near the cathode, as expected. 


Running the program over and over again with the same parameters just to get 
plots of different variables can be time consuming. So three of the plots were 
combined into one multi-plot. The multi-plot program was then modified to make it into 
a “cartoon,” that is, the plots would update on the screen with each chosen (by the 
programmer) time step. Fig. 11a shows a multi-plot cartoon freeze-frame at a time just 
before the current density reached the absolute minimum. As expected, the density 


plot is high near the anode and low near the cathode. The electric field plot shows 
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that the E-field at the half-way point across the gap is less than E,,,, also indicating 


that there truly is more charge in front of the mid-point than behind it. 


In Fig. 11b the time is frozen at a point where the current density has just 
reached a local maximum. The density graph shows high density at either end of the 
diode (newly created at the cathode and what's left of the charge at the anode from 
Fig. 11a). The electric field curve "kinks" in tandem with the new charge in the density 


plot below it. 


Fig. 11c shows the frame in which the current density approaches a local 
minimum. Note how the density is more spread out than in Fig. 11a, and how there 
is less curvature to the E-field curve. By Fig. 11d, the next local maximum of the 
current density has been reached, and by now the density is nearly equal everywhere. 
The electric field plot is near linearity. 

In Fig. 11e the multi-plot is run to the steady-state. The density is constant at all 
points except near the cathode; the electric field curve is straight and centered about 
Eapp, aS expected. Based on proper behavior of the above plots, it can be concluded 


that the simulation program is producing plausible data. 
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Figure 8a. Space charge velocity vs. distance ratio at equilibrium. 
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Figure 8b. Space charge velocity vs. distance ratio at equilibrium; first five percent of 


gap. 
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Figure 9a. Space charge density vs. distance ratio at equilibrium. 
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Figure 9b. Space charge density vs. distance ratio at equilibrium; first five percent of 


gap. 
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Figure 10a. Electric field vs. distance ratio at equilibrium. 
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Figure 10b. Electric field vs. distance ratio at equilibrium; first five percent of gap. 
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V. CONCLUSIONS 


As expected, the final equilibrium value of the current density for a pure field 
emitter was found to be less than the Child-Langmuir space charge limited current 
density for the non-run-time-constrained simulations. A comparison of the values is 
presented in Table ll. The run-time-constrained results (denoted by a ¢) cannot be 
used to either verify or contradict Eq. (4), but based on the arguments of section IID. it 


is likely that verification would be achieved if longer run-times were feasible. 


Based on the simulation results shown in Table |, it can be stated that, except in 
the cases denoted by a i, that Je<Je. Thus, for a 1 in gap, a purely field-emitted 
current density is not sufficient to cause the explosion mechanism of a whisker on the 
cathode to trigger in less than 10 ns. The above noted exceptions may or may not 
conform to the other data; computer run time constraints did not allow for a definitive 


determination. 


Again, it must be pointed out that this model of only pure field emission is quite 
elementary in that it does not include the effects of thermionic emission. A whisker 
certainly heats up due to the emission process and some thermionic emission must be 
present. In fact, based on the results of those simulations which ran until equilibrium 
was reached (MEapp S$ 3.94x10’ V/cm), it would seem that thermionic emission must 
be the dominant emission source, at least for those parameters, if the explosive model 
is truly accurate. Whether enough thermionic emission occurs to validate the 
explosive model is debatable, especially in cases such as 
m= 1, Eapp S 3.94x10° V/cm, where the field emission contributes almost nothing to 


the total emission current density and the applied electric field does not provide much 


41 


relative enhancement (compared to higher values of E,,,) at the metal-vacuum 
interface. A study done by Hallal (Ref. 3] at such parameters concluded that the 
explosive model could not account for the breakdown that he observed, and that a 
Schwirzke-type model seemed more promising. It can be concluded, then, that the 
explosive model only has a chance of validity in the case of ME,,, 2 3.94107 V/cm, 


with some small dependence on work function, for a diode with a 1in gap. 
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APPENDIX A : DIMENSIONLESS VARIABLES 


Although modern computers can easily handle both very large and very small 
numbers, it is still considered good programming practice to work with small (= 1) 
numbers whenever possible. To facilitate this, variables (such as time, distance, 
velocity, etc.) are transformed into dimensionless quantities. The transformation is 
accomplished by dividing each occurrence of a particular variable by some constant 
value of the variable relevant to the problem. For example, if one were to examine the 
motion of billiard balls on a table top, one could divide the masses in the equations of 
motion by the mass of a single billiard ball m,, thus reducing the values of the masses 
to 1 during the computations. After all of the calculations are finished, one would then 


multiply the dimensionless values by the appropriate constants to get real values. 


since some variables are composites of others (e.g., velocity is composed of 
distance and time), once a certain number of variables has been transformed, 
transformation of the remaining ones is then governed by the constants already 
chosen and the physical equations of the problem at hand. For the billiard ball 
example, one might choose time, distance, and mass for transformation first. 
Momentum, force, and energy, since they are composed of the three already 
transformed variables, would eventually be transformed using the three constants for 
time, distance, and mass, along with any physical constants which may appear in the 
governing equations. Often the "base" (non-composite) variables are chosen for 
transformation first. In this thesis, such base variables would be time, mass, charge, 
and distance. Indeed, the first three have been chosen here. But rather than 
distance, current density has been selected instead. Current density was chosen 


since the graph of current density vs. time was so central in the thesis. 


43 


The basic transformations are, then: 


time : t x 7 _) ae 
1 
m : OM — M: 
ass x Mi + 
: 1 
charge :  q x = => q 
current 
density oreere) x all ~ Aa 
Jo 


All other variables can be transformed using the four constants above: 








2M 

distance a x eS =_ x; 
if e€ Jo 
2&,M 

velocity ea’, x — 7 ==) wes 
a @ Jo 
2, M 

acceleration : a x soles = a: 
T@do 

electric 5 

field ae = ‘ — ee = 

0 
Where 


T = time step (s) ; 

M, = electron mass = 9.11x10°' kg _; 

e = lelectron charge | = 1.60x107'9C ; 

Jo = current density at cathode at t= 0 (A/cm?) ; and 


E>) = permitivity of free space = 8.85x107-'* C*/N-m? . 
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The extra factor of 2& comes from the governing physical equations. This can 


be seen by looking at the derivation of the transformation of the dimensionless electric 


field, E. Eq. (13) can be rewritten as 


Jo T Jk At = Jo T 
Hera 1 © 2, 


> Jk At. 





Rearranging gives 


2 Ep 
ail 





Ee y J, AT. 


qc 


As the right side is dimensionless, so must be the left side, and one gets 


2 & : 
E = ——E = ¥ J, At. 
ay 2 Yk 
If At = 1, as it does in the program, 
E = Y, ik : 


(23) 


(24) 


(25) 


(26) 


The power of Eq. (26) lies in its simplicity; dimensionless electric fields and current 


densities become "equivalent." Only addition is required for the calculation of the 


dimensionless electric field; no program steps are wasted multiplying 


the 


dimensionless current densities by constants. As equations of the form of Eq. (26) 


make up a large part of the computations, the reduction in the number of calculations 


can be significant. 
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Since the factor of 2& appears in the expression for the dimensionless electric 
field, it must also appear in the remaining equations for acceleration, velocity, and 


distance. From Eq. (22) one can get for the acceleration a = dv/dat 











E(2 M edo T 
goatee GENS Sol NeR Eis (27) 
yo e JoT M {|2€M, 
Rearranging yields 
2 
e Jo T ~ MA 
Again, as the right side is dimensionless, so must be the left side, and one gets 
2&M 
a = SO ec (29) 


edgT ~M 


Interestingly enough, since q = M = 1 in the program, acceleration, almost like current 


density, is nearly “equivalent” to the electric field in dimensionless form. 
Velocity and distance are related to acceleration by factors of (1/T) and (1/T)?, 
respectively: 


2 & Me 2 & Me 


V; - 
en) ealic enue 








(30) 
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APPENDIX B : SIMULATION PROGRAM 


/* Variable Parameters ‘*/ 


#include 
#include 
#include 


#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 


<stdio.h> 
<cmath.h> 
<greeks.c> 


RUN 
dt 

+ 
GAP 
m 
wkfun 
Eapp 
yscale 
MAX 
mu | t 
hunch 
limit 


100000 
1 


1. Oe 17 
2.54 
100 

4.0 
3.94E6 
1000 
ie 

‘ 

see 

800 
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* * > + FF © * > fe s 


— Ss SS TS SS SS Se SS, TS, TM 


array size */ 

time step (dimensionless) /* 
time step (sec) */ 

gap size (cm) */ 

enhancement factor */ 

work function (eV) */ 
applied E-field (V/cm) */ 

for plotting. 7 


plot all pts. w/y-val.>MAX */ 
plot limit/mult pts. */ 
artificial charge est. */ 

# of run steps */ 


/* Declarations and Assignment of Initial Values */ 


main() 

{ 

double E[RUN+1],J[RUN+1],x[RUN+1],v[RUN+1],a[RUN+1],F, 
Jmult,gapdist,dx,xmin,xmax,ymin, ymax, lablxmin, 
lablxmax,lablymin, lablymax,ploty [RUN+1],Ffact, 
Jfact,Jexpfact,c2,alpha,ganmmainv,spchrg; 

int c,f,i,j,k,1,p,q,t,oldest,n,pmult,ixmin, ixmax, 

iymin, iymax ,plotx [RUN +1],mark,iy,flag,ily; 

pmult = mult; 

Jexpfact = (-6.42E7)"pow(wkfun,1.5) ; 

Pe=-in Eappr 

J[O] = ((1.54E-6)*F*F/wkfun)*exp(Jexpfact/F) ; 


Ffract 
Jfact 


m*T*J[0]/(2.0*(8.85E-14)); 
(1.54E-6)/(wkfun*J[0]); 


E[O0] = Eapp*2.0*(8.85E-14)/(T*J[0]); 


alpha = (1.0E-6)*2.0*(8.85E-12)*(9.11E-31)/ 
(1 To (tec E-done tonne 


gapdist = alpha*GAP; 


x [0] 


O03, 
v[0] 0.0; 


c2 = alpha*Tl*(3.0E10)*alpha*T*(3.0E10) ; 
Idest = 1; 
1; 


ie 


o3-0 


=~ 


1; 
lag = 
jj = Bb 
spchrg = hunch; 
ymin = -0.001; 


d 
ark 
g OF, 
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/* Determine Electric Fields, Accelerations, Velocities, and 
Positions for Existing Disks; Eliminate Disks at Anode ‘*/ 
while(t <= limit) 
{ 
for(i = oldest; i <= (mark-1); i = i +1) 
{ 
if(i != oldest) 
Sie} = Bea) a= Sa a = SET] : 
ganmmainv = sqrt(1.0-(v[i]*v[i]/c2)); 
E[i]*“ganmainv*ganmainv*ganmainv; 


vii] + afi]*dt: 
x[i] + vfi]tdt; 


a[i] 
vii] 
x[ i] 


if(x[i] >= gapdist) 
{ 
oldest = oldest + 1; 


for(j = oldest; j <= (mark-1); j = j+1) 
E[j) = E[j] + J[oldest - 1]; 
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/* Determine Electric Fields, Velocity, Position, and 
Current Densities for Newly Created Disks at Cathode; 
Decrease Artificial Space Charge */ 


i f (flag==0) 


E{mark] = E[0O] - spchrg; 
else 

E{mark] = E[0]; 
x[mark] = x[0]; 
v{mark] = v[0O]; 


for(k = oldest; k <= (mark-1); k= k+1) 
E{mark] = E[mark] - J[k]; 


F = Ffact*E[mark]; 


J{mark] = Jfact*F"*F*exp(Jexpfact/F) ; 


spchrg = spchrg - J[1]; 


if(spchrg<0.0 && flag == 0) 
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/“ Determine Points for Plotting; Calculate Electric 
Field Changes Due to Newly Created Disk at Cathode; 


Increment Time Step */ 


if(t==pmult) 


{ 

plotx[c] = 
ploty[c] = 
Cc = c+1; 
pmult = pmult + mult; 


else if (J[mark]>MAX) 
{ 
plotx[c] 
ploty[c] 


Cc = Cc+1; 


Le 
log10(J[mark]); 


ie 
log10(J[mark]); 


if( log10(J[mark])<ymin) 
ymin = 1o0g10(J[mark] ); 


for(l = oldest; | <= (mark-1); | = 141) 
E{ I] = E{l] + J[mark]; 


t = t + dt; 


mark = mark + dt; 
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/* "Recycle" Arrays if Full; End While Loop */ 


if(mark>RUN && RUN!=! imit) 
{ 
for (q=oldest ;q<=RUN;q=q+1) 
{ 


E[q-oldest+1] = E[q]; 
J[q-oldest+1] = J[q]; 
a[q-oldest+1] = a[q]; 
v[q-oldest+1] = v[q]; 
x[q-oldest+1] = x[q]; 


} 


mark = RUN - oldest +2; 
oldest = 1; 


} 


5) 


/* Setup Graphics */ 


ymax = 0.0; 

ixmin = Q; 

ixmax = limit; 

iymin = (int)(ymin’yscale) ; 
iymax = (int)(ymax*yscale) ; 
lablxmin = 0.0; 

lablxmax = T*limit; 
lablymin = ymin; 

lablymax = 0.0; 


eae iy se ny ee ae eee ee ee ee ee ee 


pmvat( ).: 

color scale("cyanblue”") ; 

grey scale("greyscale1l"); window0() ; 
bgcol(7); erase(); color(0O) ; 


move(245,30);printf("Log Current Density Ratio vs. Time"); 
move(/75,384);printf("Log J/J");size(1,1);sub('0'); 
$ize(2,2); 

move(435,/700);printf("Time (sec)"); 


move(85,96);printf("%-3.2f", labl ymax) ; 
move(85,660);printf("%-3.2f",lablymin) ; 
move(195,682) ;printf("%3.2f",lablxmin) ; 
move(760,682) ;printf("%3.2e",lablxmax) ; 


move(810,150);printf("m = %d",m); 

move(810,250) ;greek(BETA) ; printf(" =%-2.1e",T); 

move(810,350);greek(Phi ) ; printf(" = %-2. 1f" wk fun) ; 

move (810,450) ;printf("E");size(1, 1ysubt' a’ )-sub(  p  ): 
sub(‘'p' ); s$ize(2,2); printf ("=%-2.1e",Eapp) ; 

move(810,550) ;printf(" = oe fe 

move (810,630) ;printf("old=%d" oldest): 


move(350,760);printf("J");size(1,1);sub(’0’);size(2,2); 
printf(" = %4.2e Amp-cm" ,J[0]);size(1,1); 
sup(’-’);sup('2’);size(2,2) ; 


window(224,96,800,672) ;bgcol (3) ;erase(); 
scale(ixmin, iymax, ixmax,iymin) ; 
rect(ixmin, iymax, ixmax,iymin);color(4) ; 


’ 
es oo eek ees re Ce ee ee, See ee eee * & Sea ies Ae 
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/* Plot Graph of Log Current Density vs Time; End Program */ 


catvec(ixmin,(int)(log10(J[1])*yscale) ) ; 


for(p=1;,p<c;p=p+1) 
printf("%d,%d,0,plotx[p],(int)(ploty[p]*yscale)) ; 


printt (ye 


window0();color(0);size(1,1);mode("P0O") ; 


iy = 672-(int)((672-96)*((yscale*ploty[p-1])-itymin)/ 
(iymax-iymin)); 


move(805, iy) ; 
printf("%4.3e" ,ploty[p-1]); 


ily = 672-(int)((672-96)*((yscale*ploty[1])-iymin)/ 
(iymax-iymin)); 


move(160,ily); 
printft("%4.3e" ,ploty[1]); 
} 
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